Laser-Shock Compression and Hugoniot Measurements of Liquid Hydrogen to 55 GPa 
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The principal Hugoniot for liquid hydrogen was obtained up to 55 GPa under laser-driven shock 
loading. Pressure and density of compressed hydrogen were determined by impedance-matching to 
a quartz standard. The shock temperature was independently measured from the brightness of the 
shock front. Hugoniot data of hydrogen provide a good benchmark to modern theories of condensed 
matter. The initial number density of liquid hydrogen is lower than that for liquid deuterium, and 
this results in shock compressed hydrogen having a higher compression and higher temperature than 
deuterium at the same shock pressure. 
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I. INTRODUCTION 

The properties of hydrogen at high pressure and high 
density are of great scientific interest. The equation 
of state (EOS) of hydrogen at these conditions is es- 
sential for modeling of the interior structure of gas gi- 
ant planets^"— . The large diversity in the estimation 
of Jupiter's core mass is resulted from the uncertainty 
in the EOS data especially in the region around the 
insulator-to-metal transition. The EOS of hydrogen- 
isotopes has important practical applications for iner- 
tial confinement fusion*, and metallic hydrogen is sug- 
gested as a prospective candidate of high-temperature 
superconductor^. Chemical free-energy models^"^ and 
ah initio simulation o^'^" have been used to predict the 
properties of warm dense hydrogen, but the results vary 
widely and have not converged yet. Therefore accurate 
experimental data for the hydrogen EOS are required for 
evaluation of the theoretical models and for further un- 
derstanding of the fundamental nature of hydrogen. 

It is more difficult to generate high pressures in hy- 
drogen than in deuterium because of its lower shock 
impedance. For this reason, most of the recent exper- 
imental measurements by shock compression have fo- 
cused on the heavier isotop&ii"— . However, it should 
be noted that the Hugoniots for the two isotopes do 
not scale in density. Owing to the difference in zero- 
point energy, the mole volume of liquid hydrogen is larger 
than deuterium^iii. As a result, hydrogen is expected to 
have higher compression and higher Hugoniot tempera- 
ture than deuterium at the same pressure. 

There is a large gap in the experimental achievement 
of shock compression between liquid hydrogen and deu- 



terium. The principal Hugoniot for liquid deuterium 
was measured up to 220 GPa using laser-driven shock 
waveai^. For the case of liquid hydrogen, the Hugoniot 
was studied experimentally only to 10 GPa by a gas gun 
and explosive method more than two decades agc^*'^^. 
The metallization of hydrogen on the Hugoniot is ex- 
pected to occur at much higher pressure. In this work, 
we carried out laser-shock experiments of liquid hydrogen 
to pressures exceeding 10 GPa in order to make a quanti- 
tative comparison of the hydrogen Hugoniot around the 
metal transition with the deuterium data. 



II. EXPERIMENTS 

The experiments were performed on the GEKKO 
HIPER laser facility at the Institute of Laser Engineer- 
ing, Osaka University. The laser is a neodymium-doped 
glass system operating at the third harmonics wavelength 
of 351 nm. Laser energies between 0.8 and 1.4 kJ were 
delivered to generate the shock pressures using a nomi- 
nally square pulse of 2.5 ns in duration. The laser focal 
spot of 600 /im in diameter was smoothed using Kino- 
form phase plates. This resulted in the effective laser 
intensities between 4 and 8 x 10^'^ W/cm^. 

Figure 1 shows the experimental setup and target ar- 
rangement. The cryogenic hydrogen target consisted of 
two z-cut a-quartz plates, by which the hydrogen layer 
was sandwiched. The thickness of the quartz was 50 ^m 
for both of them. These quartz plates were glued to 
kovar flanges attached to a copper cell filled with liq- 
uid molecular hydrogen at 15 K. The gap between the 
two quartz windows was 30-250 /im thick. The laser-side 
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quartz was deposited by 40 /im aluminum as a pusher. 
The aluminum layer has another important role to shield 
the quartz and hydrogen against the x-ray radiation from 
ablation plasma. The free surface of the rear-side quartz 
had an anti-reflective coating to minimize spurious re- 
flections. In this target, the quartz was designed to be 
used not only as a window material but also as a reference 
standard for the impedance-matching measurement. The 
initial density of quartz and liquid hydrogen were 2.65 
and 0.0760 g/cm'^. At the probe laser wavelength 532 
nm, the index of refraction of quartz and hydrogen were 
1.547 and 1.121—. Small change in the optical properties 
of quartz at cryogenic temperature is ignorable^^. 



III. RESULTS: SHOCK VELOCITIES 

Shock velocities in the laser-side quartz UsQ and hy- 
drogen UsK were measured using a line-imaging velocity 
interferometer system for any reflector (VISAR)^^'^^. At 
the high pressure involved in these experiments, shock 
waves in the initially transparent quartz and hydrogen 
are reflecting. Then the VISAR can provide a direct 
time-resolved measurement of shock velocities in these 
media. The probe light for VISAR was a Q-switched 
yttrium-aluminum-garnet (YAG) laser operating at a 
wavelength of 532 nm. Two VISARs were run concur- 
rently on each shot to resolve the 2tt phase-shift ambi- 
guities at shock breakout. The velocity sensitivities were 
4.14 and 14.53 km/s/fringe for quartz, and 5.71 and 20.05 
km/s/fringe for liquid hydrogen. Post processing of the 
VISAR images using Fourier transform methods deter- 
mines the fringe position to 5% of a fringe. The re- 
sulting velocities were measured to ~ 1% precision since 
shock velocities were high enough to cause multiple fringe 
shifts. 

A sample VISAR trace is shown in Fig. 2(a) and the 
accompanying profile of shock velocities in Fig. 2(b). 
The time-resolved VISAR measurements allowed veloc- 
ities to be tracked continuously during transit through 
both the standard and sample. The transit distance in 
the quartz determined from the time-integrated VISAR 
velocity agreed consistently with the thickness measured 
by interferometry before the shock experiment. Shock ve- 
locities immediately before and after the quartz-hydrogen 
interface were adopted for the impedance-matching anal- 
ysis. To precisely determine the velocity at the break- 
out, a linear fit was taken of the velocities over 0.5 ns 
and extrapolated to forward or backward the impedance- 
matching time t = 0. The velocities were calculated by 
averaging the shock profile over 0.1 ns, which is equiva- 
lent to the temporal resolution of the streak camera. 

As seen from Fig. 2(b), the shock velocity in quartz 
decreased with time more than 15% within ~ 2 ns. De- 
caying shock is a common feature of laser-driven shock 
waves. This means that the estimation of shock veloc- 
ity by transit time requires corrections taking account 
of the deceleration rate, and thus it could bring a large 



uncertainty in Us- Therefore transparent standards have 
a great advantage for laser-shock experiments, and actu- 
ally provide a significant improvement in the accuracy of 
the evaluation of the pressure P and density p along the 
Hugoniot^'^'^'*. Shock planarity was observed to extend 
over 400 /im which ensures more than 5 fringes to be 
available for the Fourier analysis. In this spatial range, 
the systematic variation in the shock arrival time was 
less than about 0.1 ns. Preheat in quartz and hydrogen 
can be neglected because preheat in the laser-side quartz 
was estimated to be less than 0.04 eV and no motion of 
the quartz-hydrogen interface was observed before shock 
breakout in ID radiation hydrodynamic simulations with 
the MULTI cod(M. 

Figure 3 shows a impedance-match diagram in the P- 
Up plane, where Up is the particle velocity. The quartz 
Hugoniot was previously established by the OMEGA- 
laser experiments to have a linear Us-Up relation given 
by UsQ = ao-t-ai(C/pQ — /3), where ag = 20.57±0.15, oi = 
1.291 ± 0.036, and /3 = 12.7423-26. xhis formula is valid 
at the pressure over 200 GPa and in good agreement with 
the Kerley-7360 EOS model^i. Impedance-matching oc- 
curs at the intersection of the hydrogen Rayleigh line 
and the release isentrope of quartz. Because of the low 
shock impedance of hydrogen, quartz release curves from 
several hundreds of GPa down to more than one order of 
magnitude are needed. In this work, the quartz isentrope 
was calculated from the tabulated Kerley model. 

Pre- and post-shock conditions are related by the 
Rankine-Hugoniot conservation equations: 

PoUs^ P (Us-Up) , (1) 
P - Po = PoUsUp , (2) 

e-eo = i(P + Po)f---) , (3) 

where e is the specific internal energy and subscript 
denotes initial conditions. The pressure and density of 
shocked hydrogen are derived from measured Us and Up 
using Eqs. ((T]) and The results for the hydrogen 

Hugoniot are listed in Table I. Random errors shown in 
the table come from measurement uncertainties in Us- 
Random error propagation in the impedance-matching 
procedure is illustrated by Fig. 3. Systematic errors aris- 
ing from uncertainties in the principal Hugoniot of quartz 
are also considered. Figure 4 depicts the principal Hugo- 
niot of liquid hydrogen in the P-p plane, where p is nor- 
malized by the initial hydrogen density pQ. These results 
are shown along with the past experimental datapi^iiS and 
theoretical models calculated for hydrogen with the ini- 
tial density po = 0.071 g/cm^^J^. 

The pressure range of shocked hydrogen obtained in 
this work was 25-55 GPa. The highest pressure is more 
than 5 times higher than the previous data. The com- 
pression p/ Po of hydrogen ranges from 3.8 to 5.6, which 
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is mostly comparable to the models but slightly softer at 
30-50 GPa. Note that the slope change near 15 GPa in 
the Hugoniot curve predicted by the quantum molecular 
dynamics simulations is related to the dissociation of hy- 
drogen molecules^^. The compression of deuterium, on 
the other hand, has been measured as 3.3-4.4 at the same 
pressure rangei^ii^"^, so that hydrogen exhibits system- 
atically higher compression compared to deuterium (see 
Fig. 4). This trend is qualitatively consistent with the- 
oretical predictions^ and can be seen in the earlier data 
below 10 GPaiS. 

The uncertainties in the Hugoniot pressure and density 
are caused mostly by the random errors associated with 
velocity measurements, 6Us and Slip, which are evaluated 
by 




As seen from Fig. 4, the uncertainties in density are 
much larger than those in pressure. Since hydrogen is 
highly compressible at this pressure range, the factor in 
Eq. dS]) becomes considerably large, p/po — 1 ~ 4. Then 
the errors in density are enhanced up to about 6%, even 
when the velocity errors are sufficiently small, SUs/Ug ^ 
SUp/Up - 1%. 

IV. RESULTS: TEMPERATURE 

Using a streaked optical pyrometer (SOP), the 
temperature T of shocked hydrogen was measured 
simultaneously^ii^. Despite temperature is fundamen- 
tal to thermodynamics, T is not a part of the Rankine- 
Hugoniot relations, and thus must be measured sepa- 
rately from shock pressure and density. We extracted 
the temperature by fitting the absolute spectral radiance 
I\ to a gray body Plank spectrum, 

/ , X 2Trhc^ 1 , , 

h{T)^e(X) (6) 

with T as a fit parameter, where c is the speed of light, 
h is Planck's constant, and k is Boltzmann's constant. 

The radiance was obtained in a single spectral band 
centered at 455 nm with a bandwidth 38 nm. The emis- 
sivity £ is related to the reflectivity R of the shock front 
through e(A) = 1 - R{\). For the reflectivity at the SOP 
wavelength, we adopted the reflectivity measured at the 
VISAR wavelength assuming a weak dependence of R 
on the wavelength (see Fig. 7). To relate the diagnos- 
tic system output to a source radiance, in situ measure- 
ments were performed to determine the spectral response 
of the system. We used calibration signals recorded with 



a 3000-K tungsten quartz-halogen reflectorized lamp to 
find the radiance at our SOP wavelength. The radiance 
of this lamp was calibrated with a bolometric calibra- 
tion traceable to the National Institute of Standards and 
Technology (NIST). 

Absolute reflectivity is determined by comparing the 
shock reflectivity to that from the aluminum surface 
which has a known value, 86 ± 6%. The obtained R of 
shock front in hydrogen is listed in Table II. The reflectiv- 
ity just after the shock breakout was calculated by the 
same algorithm as the shock velocity Ug using VISAR 
signals. The reflectivity of shocked hydrogen increases 
from 8% to 23% as the pressure increases from 25 GPa 
to 41 GPa. The Fresnel formula for the reflectivity of 
the shock front is i? = \{hs — no)/(ns -I- jt-o)!^: where hg 
is the complex refractive index behind the shock front 
and no is the refractive index in the undisturbed liquid. 
At low pressures, the refractive index of liquid hydrogen 
follows an empirical form, n — 1 cx p^. Assuming the 
fluid remains mostly in its molecular form under shock 
compression, the refractive index at fivefold compression 
would be n ~ 1.6, and the reflectivity ^ 4%. The much 
higher reflectivity we observed suggests that the fluid be- 
comes conducting. 

The Hugoniot temperature can be obtained using the 
reflectivity and SOP spectral intensities. An example of 
the intensity versus time data is shown in Fig. 5. The 
emission dropped dramatically when the shock front en- 
tered the hydrogen layer from quartz. Compared to the 
VISAR data, material boundary on SOP records is not 
clear due to the poorer temporal resolution ^ 0.2 ns. 
The self-emission profile of hydrogen just after the shock 
breakout is then determined from extrapolating a linear 
fit to measurements back to t ~ 0, in order to eliminate 
the contamination of the quartz emission. The time aver- 
age of the fitted intensity over the temporal resolution is 
adopted for the derivation of the Hugoniot temperature 
of hydrogen. The results of the hydrogen temperature are 
listed in Table II. The size of the uncertainty in temper- 
ature is 6-19%, which stems from the system calibration, 
optical transmission, and measurements of self-emission 
and reflectivity. 

Figure 6 depicts the Hugoniot temperature of liquid 
hydrogen as a function of measured pressure. Theoreti- 
cal models predicted that Hugoniot temperatures of liq- 
uid hydrogen are higher by a factor of 1.3 than those of 
liquid deuterium at the same pressure- . Deuterium tem- 
perature measured by previous shock experiments^!^^ 
are also shown in Fig. 6. Although the uncertainty 
in temperature measurements is still large, the hydro- 
gen temperature obtained in this work is higher than the 
deuterium data. This is another experimental evidence 
of the isotope difference in the hydrogen Hugoniot. At 
this temperature, molecular hydrogen begins to dissoci- 
ate and the fraction reaches to ~ 14% (39%) at P = 25 
GPa (40 GPa) according to Kerley's model^. 
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V. DISCUSSION AND CONCLUSIONS 

Drude-type models are often applied to parameterize 
the optical properties of liquid metal^^i^. Within the 
Drude description, the complex index of refraction is 
given by = 1 — (Wp/w^)(l + i/ojTf.)~^ where ujp = 
(47rnee^/me)^/^ is the plasma frequency, rie is the carrier 
density, e is the electron charge, mg is the electron mass, 
and u) — 2-Kc/ A is the optical frequency. Here the electron 
relaxation time is assumed to be = Rq/vf where Rq 
is the interparticle spacing and is the electron Fermi 
velocity. We adopt Rq — 0.126 nm corresponding to a 
density 0.4 g/cm'^. Then the Drude reflectivity of liq- 
uid metallic hydrogen can be derived as a function of rie, 
which is depicted in Fig. 7. 

High reflectivity is produced when the carrier density 
exceeds the critical density ric defined by ujp{nc) = uP' . 
The critical density is Uc = 3.9 x 10^^ cm^'^ at A = 532 
nm. Taking R = 23%, for example, the carrier density 
is given by 3.7 x 10^^ cm"'^, which is about 31% of total 
hydrogen number density rimax = 1.2 x 10^"^ cm^'^ and 
equivalent with = 0.37. The Fermi energy estimated 
from this carrier density is 4.0 eV, which is higher than 
the shock temperature. This indicates that the observed 
highly reflective state would be characteristic of a degen- 
erate liquid metal. Based on the Drude model, the reflec- 
tivity increases with the carrier density and saturates at 
R = 36%. This implies that the maximum reflectivity of 
hydrogen may be slightly lower than that of deuteriumSi^ 
although it depends on the assumed compression ratio. 
It will be interesting to confirm the saturated value of 
the reflectivity by future experiments. 

The hydrogen Hugoniot obtained in this work is 
strongly dependent on the quartz EOS. Recently the 
quartz Hugoniot was reexamined by using magnetically 
driven flyer impact on the Z machine and a new fitting 
function was derived'^^. For comparison, the hydrogen 
Hugoniots calculated by using the Z-fit and Kerley re- 
lease (J7pjj, Py, and p^) are listed in Table I. The stiffer 
Hugoniot of quartz reduces the initial pressure of the re- 
lease isentrope (see Fig. 2). At the initial shock state of 
quartz, the pressure and density inferred from the Z-fit 
Hugoniot, Pi — 716 GPa and pi — 6.61 g/cm'^, are lower 
than those of the OMEGA-fit case. Pi = 746 GPa and 
pi = 7.05 g/cm^). 

In our analysis, the quartz release cannot be approxi- 
mated by a reflection of the Hugoniot in the P-U^ plane, 
and thus the off-Hugoniot EOS of Kerley's model deter- 
mines the shape of the curve. The particle velocity along 
the release isentrope can be calculated by using the Rie- 
mann integral: 

Up = C/pi - I Csp~^dp , (7) 

where Vp\ is the particle velocity at the initial shock 
state, and the sound speed is obtained from the pres- 
sure derivative of the density at constant entropy S", 



= {dP/dp)s- The second term of Eq. ([7]) is larger 
for the Z-fit case due to the lower density, and then the 
particle velocity increases faster with the decrease of den- 
sity. Therefore, two release curves shown in Fig. 2 are 
gradually approaching at the lower pressure. The result- 
ing differences in the Hugoniot pressure and density for 
hydrogen are 1% and ^ 4% compared to those de- 
rived by the OMEGA-fit. Apparently, improvement of 
the quartz EOS is essential for the further development 
of Hugoniot measurements using quartz standards^S. 

In summary, we have obtained the principal Hugoniot 
data P-p-T for liquid hydrogen, not deuterium, in an 
unexplored range of pressure up to 55 GPa. The re- 
sults demonstrate that the hydrogen Hugoniot cannot 
be scaled by density from the deuterium data in conse- 
quence of the initial density effectsi^i^li^S. As for the 
study of planetary interiors, the hydrogen EOS data at 
much higher pressure are required since the transition to 
metallic hydrogen is anticipated to be at P 200-400 
GPa in Jupiter™. However, the hydrogen temperature 
must be kept lower because the Hugoniot temperature 
at this pressure range is too high to reproduce Jupiter's 
conditions. Therefore off-Hugoniot measurements of hy- 
drogen by means of reflection shocks^ and/or precom- 
pressed samples^"— will be a quite important next step. 
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FIG. 1: (Color online) Sketch of the cryogenic hydrogen target 
used in the experiments. The drive laser irradiates the target 
from the left with an incident angle of 30 degree, while the 
VISAR and SOP measure the shock velocity and self-emission 
from the rear-side of the target. 
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FIG. 3: (Color) Impedance-matching in the P-Up plane il- 
lustrating measurement error propagation. Black dot-dashed 
curves show the Rayleigh lines for quartz (upper) and hydro- 
gen (lower). Red solid curve is the quartz Hugoniot derived 
by laser experiments (Ref. 26). Blue solid curve is a release 
isentrope of quartz calculated from Kerley's model (SESAME 
7360) (Ref. 27), and blue dashed curves show the propaga- 
tion of uncertainty in the quartz shock velocity. The inset is 
a magnified view around the intersection of the Rayleigh line 
of hydrogen with the release isentrope of quartz. Black dot- 
ted curves show the range of uncertainty in the Rayleigh lines 
associated with the measurement error in hydrogen shock ve- 
locity. For the purpose of comparison, the quartz Hugoniot 
based on Z experiments (Ref. 35) and the corresponding Ker- 
ley release are shown by orange and green solid curves. 
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FIG. 2: (a) Sample VISAR trace showing the signal from the 
reflecting shock front in the laser-side quartz and liquid hydro- 
gen in experiment No. 32275. (b) Resulting velocity profile 
extracted from the VISAR fringe shift in (a) with dashed lines 
representing measurement uncertainties. 



TABLE I: Hugoniot data for liquid hydrogen resulted from impedance-rnatcliing to a quartz standard. UsQ and Usu are the shock Aclocities in quartz and hydrogen 
with measurement errors. Upu, -Ph, and pu are the particle velocity, pressure, and density of shocked hydrogen showing both random and systematic errors: (ran, sys). 
Upii, Ph, and pn are reference data derived by using the Z-fit quartz Hugoniot (Ref. 35) instead of the OMEGA laser fit (Ref. 26). Shot numbers with an asterisk 
indicate targets with a thinner aluminum pusher (23 fj,m) and VISARs with lower sensitivities, 9.79 (13.51) and 14.72 (20.31) km/s/fringe for quartz (hydrogen). 



Shot No. 


UsQ [km/s] 


UsH [km/s] 


UpH [km/s] 


Ph 


[GPa] 




PH [g/cm^] 




Ku [km/s] 




P' 


[GPa] 


Ph [g/cm='] 


32275 


21.24 ±0.21 


30.31 ± 0.20 


24.16 ±(0.38, 0.01) 


55.7 ± 


(0.9,0. 


,0) 


0.375 ± (0.024, 0.001) 


24, 


.03 ± (0.39, 0. 


,14) 


55.4 ± 


(1.0,0.3) 


0.367 ± (0.023, 0.008) 


33159 


19.33 ±0.16 


26.89 ± 0.30 


21.44 ± (0.27,0.04) 


43.8 ± 


(0.7,0. 


■1) 


0.375 ± (0.025, 0.003) 


21 


.29 ± (0.29, 0. 


,13) 


43.5 ± 


(0.7,0.3) 


0.365 ± (0.025, 0.008) 


33144 


19.27 ±0.14 


26.64 ±0.20 


21.37 ± (0.26,0.05) 


43.3 ± 


(0.6,0, 


.1) 


0.384 ± (0.020, 0.003) 


21 


.21 ± (0.27,0, 


,13) 


42.9 ± 


(0.6,0.3) 


0.373 ± (0.020, 0.009) 


32270 


19.04 ± 0.34 


25.71 ± 0.25 


21.10 ± (0.54,0.06) 


41.2 ± 


(1.1,0, 


.1) 


0.424 ± (0.054,0.005) 


20 


.92 ± (0.56, 0, 


,13) 


40.9 ± 


(1.2,0.3) 


0.408 ± (0.054,0.011) 


31912* 


17.54 ±0.28 


23.72 ± 0.45 


18.95 ±(0.49, 0.07) 


34.2 ± 


(1.1,0, 


■1) 


0.378 ± (0.053, 0.006) 


18 


.77 ± (0.44, 0. 


,20) 


33.8 ± 


(1.0,0.4) 


0.364 ± (0.045, 0.014) 


33151 


16.58 ± 0.17 


22.21 ±0.27 


17.43 ±(0.32, 0.11) 


29.4 ± 


(0.6,0, 


.2) 


0.353 ± (0.028, 0.008) 


17 


.50 ±(0.28,0. 


,22) 


29.5 ± 


(0.6,0.4) 


0.358 ± (0.027, 0.016) 


31917* 


15.49 ± 0.44 


21.14 ±0.54 


15.78 ±(0.68, 0.14) 


25.4 ± 


(1.2,0, 


■2) 


0.300 ± (0.052, 0.008) 


15 


.97 ± (0.67,0, 


,17) 


25.7 ± 


(1.2,0.3) 


0.311 ± (0.054,0.010) 


31922* 


15.26 ± 0.48 


20.93 ± 0.71 


15.45 ±(0.74, 0.14) 


24.6 ± 


(1.4,0, 


■2) 


0.290 ± (0.062, 0.008) 


15 


.64 ± (0.74, 0, 


,14) 


24.9 ± 


(1.4,0.2) 


0.301 ± (0.066, 0.008) 
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FIG. 4: (Color online) Pressure versus compression for the 
principal Hugoniot of liquid hydrogen. Data are from Dick 
and Kerley (Ref. 18) (blue triangles), Nellis et al. (Ref. 19) 
(green squares), and this work (red circles). Error bars rep- 
resent the quadrature sum of random and systematic errors. 
Also shown are theoretical predictions from the model EOS 
of Kerley (Ref. 8) (solid line), quantum molecular dynamics 
simulations (Ref. 10) (dashed line), and the linear mixing 
model (Ref. 28) (dot-dashed line). For reference, the Hugo- 
niot data for liquid deuterium (Refs. 14, 16, and 30) are de- 
picted by gray diamonds, and gray curve is the Kerley model 
for deuterium (Ref. 8). 
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FIG. 5: Time profile of the emission intensity at 455 nm cal- 
culated from the central region of the streaked shock emission 
data for experiment No. 33151, which is shown in the inset 
image. The thickness of the hydrogen layer was ~ 250 /^m for 
this shot, so that the transit time is much longer than 2.5 ns. 
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FIG. 6: (Color online) Temperature versus pressure for the 
principal Hugoniot of liquid hydrogen. Experimental data are 
from Holmes et al. (Ref. 30) (green squares) and this work 
(red circles). For comparison, the Hugoniot temperature for 
liquid deuterium (Refs. 30 and 32) are depicted by gray dia- 
monds. Also shown are theoretical predictions from Kerley's 
EOS model (Ref. 8) for hydrogen (solid line) and deuterium 
(gray line). 
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TABLE II: Reflectivity iin and Hugoniot temperature Th 
for liquid hydrogen. The reflectivity of shocked hydrogen was 
measured at the wavelength 532 nm using VISAR signal and 
aluminum known reflectivity. The Hugoniot pressure and den- 
sity for each shot are listed in Tabic I. 



Shot No. 
32270 
33151 
31922 



i?H [%] 

22.7 ±2.8 
14.1 ± 1.9 
8.3 ±1.8 



Th [eV] 
0.93 ±0.08 
0.77 ±0.05 
0.69 ±0.13 



FIG. 7: Drude reflectivity plotted as a function of carrier 
density at the VISAR wavelength 532 nm (solid line) and 
SOP wavelength 455 nm (dashed line). The difference in the 
reflectivity between these two wavelengths is at most 2%. The 
critical density ric and maximum density Wmax of electrons are 
indicated by the arrows. 
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